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Dynamical evolution of test fields in background geometry with a naked singularity is an important 
problem relevant to the Cauchy horizon instability and the observational signatures different from 
black hole formation. In this paper we study electromagnetic perturbations generated by a given 
current distribution in collapsing matter under a spherically symmetric self-similar background. 
Using the Green’s function method, we construct the formula to evaluate the outgoing energy flux 
observed at the future null infinity. The contributions from “quasi-normal” modes of the self-similar 
system as well as “high-frequency” waves are clarified. We find a characteristic power-law time 
evolution of the outgoing energy flux which appears just before naked singularity formation, and 
give the criteria as to whether or not the outgoing energy flux diverges at the future Cauchy horizon. 
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I. INTRODUCTION 

While Penrose has proposed the so-called cosmic censorship conjecture^l, several spacetimes admitting naked 
singularity formation even in gravitational collapse of physically reasonable matter from regular initial data have been 
found over the past three decades (see for a recent review). The examples well-studied are spherically symmetric 
collapses of an inhomogeneous dust ballP, 0, and a self-similar isothermal gas|^ S IE 0' ^ shell-focusing naked 
singularity can appear at the center in a wide range of the initial data set. 

It is an important problem to understand peculiar physical phenomena associated with naked singularity forma¬ 
tion. Time evolution of various perturbations in black hole geometry have been extensively studied, and the late-time 
behaviors such as quasi-normal ringings and power-law tails have been clarified. Such perturbation analysis in back¬ 
ground geometry involving a naked singularity will be interesting not only in relation to the possible instability of 
the Cauchy horizon but also in terms of revealing some typical patterns of time evolution observable as a precursor 
of naked singularity formation. 

The first step to approach this problem will be to consider spherically symmetric self-similar models of gravitational 
collapse as background spacetime. Because all dimensionless components of the background metric depend only on 
the one self-similar variable z = r/t, the perturbation analysis becomes mathematically simpler. In addition, there 
exist numerical simulations ^3 showing that the geometrical structure and the fluid motion at late stages in general 
non-self-similar collapse of an isothermal gas can be well described by the general relativistic versionj^ of the Larson- 
Penston self-similar solution. We can expect self-similar solutions to be a realistic model of gravitational collapse. 

Several works have been devoted to the analysis of perturbations generated in spherically symmetric self-similar 
collapse with naked singularity formation. For example, using quantum theory of particle production in curved 
spacetime, it has been shown that the energy flux of the semiclassical radiation diverges on the Cauchy horizon 
according to inverse square power-law of the retarded timeflll Ha Hq . On the other hand, considering 

classical perturbations of a massless scalar field, Nolan and WatersflTj have claimed that the energy flux of scalar 
perturbations remains finite even at the Cauchy horizon. Their analysis is based on the behavior of perturbations 
written by v'^Hn{z) with the advanced null coordinate v and a complex parameter n. In this paper we pursue 
the classical analysis more extensively, by using the Green’s function method. Though we use a standard Fourier 
decomposition of the Green’s function by the function exp(ia; log |t|) (in an analogous way to the Mellin transformation 
in Ea ), we reveal time evolution of perturbations by integrating the Fourier components with respect to the spectral 
parameter w, which will allow us to find new features missed in EE 

Further we focus our investigation on electromagnetic perturbations generated by any given axisymmetric and cir¬ 
cular current distribution in collapsing matter. Of course we must assume that the existence of the source current 
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does not disturb the background self-similar metric. This special choice of a test field is partly motivated by an astro- 
physical interest related to highly energetic phenomena such as 7 -ray bursts (for example, see la for electromagnetic 
radiation in a process of black hole formation). In addition, the mathematical setup of the Green’s function method for 
electromagnetic perturbations becomes quite simpler in comparison with gravitational perturbations. Nevertheless, 
we would like to mention that the Green’s function method developed in the following sections does not essentially 
rely on a special property of electromagnetic fields. The application to scalar and gravitational perturbations will be 
straightforward. 

In this paper we are mainly concerned with time variation of the outgoing energy flux (namely, the Poynting flux) 
observed at the future null infinity. In Sec. nn we illustrate spherically symmetric self-similar models admitting naked 
singularity formation at the center, and we derive the basic equation for a gauge invariant and odd parity electromag¬ 
netic perturbation. In Sec. cni we introduce the retarded Green’s function and its Fourier decomposition, from which 
we construct the formula to extract the outgoing wave part of the electromagnetic perturbation propagating to the 
future null infinity. Contributions from various Fourier components parameterized by lo, (for example, corresponding 
to high-oj waves, quasi-normal modes of the self-similar system with a complex uj) are clarified in Sec. HU A charac¬ 
teristic power-law behavior (accompanied with or without an oscillatory behavior) of the outgoing energy flux as a 
function of the retarded time is found as a signature of naked singularity formation. We obtain in the final section 
the critical conditions for the self-similar background as to whether or not the outgoing energy flux diverges at the 
future Cauchy horizon, and a physical interpretation of our results is presented. Throughout this paper, the units in 
which c = G = 1 are used. 


II. SETUP OF THE SYSTEM 


Let us begin with a brief review of spherically symmetric self-similar collapse with naked singularity formation. 
The line element of the background spacetime is given by 

ds^ = -h r^S^{d9^ -b sin^ Odcj)^), (2.1) 

with the comoving coordinates t and r. The self-similarity considered here means 

v = v{z), A = A(z), S = S{z), (2.2) 

where z = rlt, and the Ricci scalar has the form 


R = — X (function depending only on z), 


(2.3) 


which should grow up as t approaches zero along a, z = constant line. Thus a singularity will appear at the center 
r = 0 at the time t = 0. A spacelike singularity may also appear along the z = z^, line where 5'(z*) = 0. 

Using the function defined by 


V{z) = ze^-', 


(2.4) 


the equations for the radial outgoing and ingoing null geodesics are given by 


dz z{l ± V) 

Hi ^ tv 


(2.5) 


The key point is that the z = constant line giving V = 1 (or V = —1) becomes a radial outgoing (or ingoing) null 
geodesic. In order to allow the central singularity to be naked as a result of collapse with regular initial data, in 
this paper, we assume the schematic behavior of the function V(z) drawn in Fig.^ for which we obtain the global 
structure of the self-similar background written in Fig. |21 (see Q for details of how to construct kinematically the 
spacetime diagram). 

The value of local minimum of U at z = Zm should be less than unity, because the line z = Z 2 > Zm giving V = 1 
becomes the future Cauchy horizon for the null naked singularity corresponding to another root z = Z 3 < Zm of V = 1. 
There exist the future null infinity and the event horizon also at the z = Z 2 and z = Z 3 null lines, respectively. The 
equation U = — 1 has a unique root z = zi < 0 corresponding to the past Cauchy horizon and the past null infinity. 
Such a behavior of U as a function of z is found in a wide parameter range of the self-similar metrics describing the 
collapse of dust or isothermal gas miiiMiii. 



3 



FIG. 1: Schematic description of the function V(z) giving a typical collapse with naked singnlarity formation. In this figure it 
should be noted that y(0) =0 and V —> ±oo as a —> ±oo. 


singularity 

Z= Zf 



FIG. 2: The Penrose diagram of spacetime corresponding to the function V(z) shown in Fig.Q The values of the null coordinates 
u and V are given at some boundaries of the regions I, II and III. 


We are interested in time evolution of electromagnetic perturbations in the regions I, II and III shown in Fig. 
For this purpose it is useful to introduce the retarded and advanced null coordinates u and v. Integrating Eq. (12.511 . 
we obtain 

U =(2.6) 

where u < 0 in the region I and n > 0 in the region II and III, while u < 0 in these regions. The functions /lout(-z) 
and hin{z) are expressed by 


f Vdz f Vdz 

hout{z)= j h,^{z) = j (2.7) 

with the requirement that /lout = Kn = 0 at z = 0 (namely, u = v = t<0a,t the regular center r = 0 represented 
by the z = 0 line). We can also note that the null coordinates are given by m = u = r at the spacelike hypersurface 
t = 0 corresponding to the horizontal z = ±oo line, which is the boundary between the regions II and III. In the limit 
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z ^ zi the function h-m is approximately given by 

--^^log|z-Zi|. (2.8) 

ziV'{zi) 

Hereafter a prime means the derivative with respect to z. Because —zit^'(zi) is positive (see Fig.^l, the advanced null 
coordinate v becomes continuously zero, even if not analytically, at the boundary z = Zi (namely, the past Cauchy 
horizon) between the regions I and 11. It remains finite at the past null infinity where z = zi and t u = —oo. On 
the other hand, in the limit z Z 2 'Vfe obtain 

^out(2) ^^^;7^^1og|z-zal- (2.9) 

Because Z 2 V'(z 2 ) is also positive, the retarded null coordinate u also becomes continuously zero at the boundary 
z = Z 2 (namely, the future Cauchy horizon) of the regions III, though it remains finite at the future null infinity where 
z = Z 2 and t ~ u = oo. It is easy to check that in the limit ^ = A = 0, the coordinates u and v become identical with 
the usual null coordinates giving the Minkowski metric. 

The values of —ziV'(zi) and Z 2 V'(z 2 ) in Eqs. (12.811 and (12.911 will become important to determine the behavior of 
the outgoing energy flux of the perturbations near the future Cauchy horizon. In the spherically symmetric spacetime 
written by the self-similar metric (lOl . there exists a homothetic Killing vector field = (—t, —r, 0,0), which satisfies 


^x9ab = -2gab, ( 2 . 10 ) 

where denotes the Lie derivative along the field x“ (the minus sign of the t and r components of x“ is chosen so 
that x“ points toward the direction toward which the curvature increases^). The values of —ziV'(zi) and Z 2 V'{z 2 ) 
are derived by evaluating the derivative of the inner product of the homothetic Killing vector x“ on the past Cauchy 
horizon z = zi and the future Cauchy horizon z = Z 2 as follows, 

V“(x'a6) = -2kx“, (2.11) 


where k = —ziV'{zi) on the past Cauchy horizon and k = Z 2 V'{Z 2 ) on the future Cauchy horizon. In the stationary 
black hole geometry, the constant k evaluated on the event horizon for a Killing field x“ normal to the horizon has 
a physical meaning as the surface gravity. Though —ziK'(zi) and Z 2 V'{z 2 ) do not necessarily have such a physical 
meaning, in this paper we denote — ziK'(zi) and Z 2 K'(z 2 ) by ki and K 2 , respectively. 

Finally in this section we derive the basic equation for electromagnetic perturbations in the self-similar background, 
by considering the Maxwell equation \7^{dbAa — da At,) = 4:Trja for the vector potential Aa- The given source current 
ja satisfying the continuity equation = 0 is assumed to be axisymmetric and circular, and in order to separate 

the variable 6 , we resolve the source current into spherical harmonics. Then, the multipole component parameterized 
by ^ = 1 , 2 , • • • is given by 


ja = (^ 0 , 0 , 0 , ji{t,r) sinO^'^ , 

which generates the vector potential Aa (with the angular momentum parameterized by ^) of the form 

Aa = ^ 0 , 0 , 0 , a(<,r) sin , 

according to the equation 

La{t,z) = -Qi{t,z), 


where 


f 2 9^ r. 9^ 


- e 


2{u-\) 




— z 


dz 


, / a €^•' 1(1 + 1 ) 

-(A + 


( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


and 


Qi{t,z) = 


( 2 . 16 ) 
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The current distribution ji(t, r) in Eq. 12.1211 is not specified in this paper, except that it rapidly decreases to zero at 
the regular center r = 0, at spatial infinity r —>■ oo and at distant past t —>■ —oo. We will obtain the solution a(t, z) 
under the boundary condition such that a(t, z) remains regular at the regular center z = 0 (and t < 0), and there is no 
ingoing flux from the past null infinity z = Zi (and t —> —oo). The latter is required to assure us that electromagnetic 
perturbations are generated through the process of the self-similar gravitational collapse. 


III. THE GREEN’S FUNCTION METHOD 


As mentioned in Sec. Q] the main purpose of this paper is to estimate the outgoing energy flux (the Poynting 
flux) observed at the future null infinity near the future Cauchy horizon, by solving the wave equation 12.1411 for the 
electromagnetic vector potential a{t,z). Our approach to this problem is to develop the Green’s function method 
useful for the analysis of any test fields in the self-similar background. 

Following the standard technique, we give the time evolution of a{t, z) as 

a{t,z) = J J G{t,z\s,y)Qi{s,y)dsdy, (3.1) 

with the (time-domain) Green’s function G{t, z\s,y), satisfying 

LG{t,z\s,y) =-d{t - s)5{z - y). (3.2) 


The causality condition requires that G = 0 if the source point represented by the coordinates s and y is located in 
the exterior of the past light cone of the observation point represented by t and z, in other words, ii u' > u or v' > v 
for the null coordinates defined by u' = u(s,y) and v' = v{s,y). 

The key point for analyzing the Green’s function is that the differential operator L in Eq. is invariant under 
the scale transformation t kt with a parameter fc, for which z is fixed. This motivates us to adopt the standard 
Fourier decomposition of the form 

Git,z\s,y) = —J G(z,y,a;)G“i°8l‘/^lda;, (3.3) 

because Eq. (Iin|l reduces to the tractable ordinary differential equation 


G"iz,y,uj) +p{z,ui)G'{z,y,ui) + q{z,ui)G{z,y,uj) 


6{z-y) 


(3.4) 


where the coefficients are given by 

, , 2(iw-l)z-f (^'-A')e2G-A)(i_^y2) 

- V^) -' 

and 

ui{i + uj)S'^z'^ — iuj{v — \)S^ z^ — + 1) 

'?G,W) = g2^2^2{.-X)^l_Y2) 


(3.5) 


(3.6) 


The spectral parameter oj introduced here may be regarded as a wave frequency in the lapse of the logarithmic time 
log|t|. Hence, hereafter we will call G{z,y,uj) the frequency-domain Green’s function as usual. 

It is straightforward to construct the frequency-domain Green’s function by the help of two independent homoge¬ 
neous solutions for Eq. ea. If the boundary condition mentioned in the previous section is imposed on a(t, z), one 
of the homogeneous solution denoted by iIjo{z,uj) should be regular at the regular center z = 0“, while the other 
denoted by ipo^t^z^uj) should be purely outgoing at z = Zi < 0 to assure the absence of ingoing waves originated from 
the past null infinity. Then, we have 


f Gi{z,y,uj) = ijjoutiz,u;)'tpo{y,uj)/W{y,uj) for z < y, 
I G 2 {z,y,uj) = ilJo{z,uj)ipout{y,u;)/W{y,u}) for z > y, 


(3.7) 


where the common factor W{y, lo) 

V’outV’o 


= w{uj)eGy)->'G )jg derived from the Wronskian 
^utV’o = w(w)exp j p(z',w)dz'^ = g 2 (.-A)^ _ v2) ’ 


(3.8) 
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with lu(uj) independent of z. 

If a homogeneous solution which becomes purely ingoing at z = Zi is denoted by '0in(z,iz>), the mode ipQ regular at 
z = 0 “ may be written by the sum 


= '0out(2:,w) +'!/'in(2:,w). (3.9) 

In this paper we claim the wave modes ij^out and tpin to become purely outgoing and ingoing, respectively, at z = zi 
in the sense that they can be expressed by the WKB forms 


V'out(-z,w) = gout{z,uj)exp{iujhoutiz)} , '0in(z,w) = gin{z, oj) exp {iujhin{z)} , (3.10) 

with the amplitudes gout and g-m “regular” at z = zi (where V = —1 and hm —oo). However, in general, we cannot 
expect that gout and gin remain regular at another regular singular point z = Z2 (where V = I and hout —^ —oo) in 
Eq. 113.411 . For example, gout may contain a term with the oscillatory factor exp(—itu/iout) which becomes singular in 
the limit z ^ Z 2 . This corresponds to a partial conversion of outgoing waves into ingoing ones in the propagation 
from the past Cauchy horizon z = zi toward the future null infinity z = Z2, and may be interpreted as a result of 
back scattering effect of spacetime curvature. (In the later discussion we will take account of this mode conversion 
to estimate electromagnetic radiation observed at the future null infinity.) Then, denoting two independent modes 
which become purely outgoing and ingoing at z = Z2 by ' 0 out and 'ipin, respectively, we obtain 

tpoutiz,u;) = bi{uj)'tpoutiz,uj) + b2{uj)'4)iniz,uj), ( 3 . 11 ) 

with some coefhcients bi and 62 dependent on to. These new modes V’out and ipin also can be written as 

V'out = gout( 2 :, w) exp{*a;/iout(^)} , V>in = gin(- 2 , tu) exp {iw/ii„(z)} , ( 3 . 12 ) 

with the amplitudes gout and gin regular at z = Z 2 , and the relation 

i’o(z,Uj) = ^ont{z,io) +'lj}in{z,Uj), ( 3 . 13 ) 


is also assumed for their normalization. 

Now let us present a more definite expression of the double integral in Eq. EU over the regions I, II and III shown 
in Fig. [5] Recall that the frequency-domain Green’s function G(z, g, to) is constructed under the conditions at the 
inner boundary z = 0 and the outer one z = zi surrounding the region I. Then, if the observation point (t, z) is 
located in the region I, the solution a(t, z) can be written as 

nZ nSx pZi nS2 

a{t,z)= dy ds Gi{t,z\s,y)Qi{s,y) + dy ds G 2 {t, z\s,y)Qi{s,y), (3.14) 

^0 J—00 Jz — 00 

where Gi{t, z\s,y) and G 2 (t, z\s,y) are given by the Fourier-type integral 13.311 of the frequency-domain Green’s 
functions Gi(z,y,uj) and G 2 {z,y,uj), respectively. From the causality condition the upper limits si{u,y) and S 2 {v,y) 
in the integrals with respect to s should be determined by the relations u(t, z) = u'{si,y) and v(t, z) = v'{s 2 , y)- Note 
that the integral of the second term of the right hand side with respect to s approaches the integral in the range from 
s = —00 to s = t as z ^ Zi (and u —> 0 ). 

The integral representation iimii assures the absence of unphysical ingoing radiation which may make the solution 
a(t,z) singular at the past Cauchy horizon z = zi (or v = 0), where the second term (including G 2 ) vanishes, and 
owing to the relation Gi ~ ' 0 out( 2 :,w) the first term (including Gi) becomes regular. This regularity of a(t,z) at the 
outer boundary z = zi of the region I allows us to extrapolate Eq. (imi to the form valid at the observation point 
{t, z) located in the region II where —00 < z < Zi and f < 0 as follows, 

pZi pS\ 

a{t,z) = dy ds Gi{t,z\s,y)Qi{s,y) 

Jo J — 00 

+ f dy f ds {Gi{t,z\s,y) - G 2 it,z\s,y)}Qi{s,y). (3.15) 

J Zi J S2 

The first term in Eq. 13.1511 corresponds to the contribution from the source point (s, y) located in the region I, while 
the second one corresponds to the contribution from the source point located in the region II. Note that the integral 
of the second term with respect to s is limited to the range S 2 < s < Si by virtue of the causality condition, and 
approaches the integral in the range from s = —ootos = tasz^zi (and u —> 0 ) to allow us to continuously shift 
from the integral representation (imi to that of Eq. 13.1511 beyond u = 0 . 
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The next step is to extrapolate Eq. to the form valid in the observation point {t, z) located in the region III 

where +oo > z > Z 2 and t > 0. As shown in Fig. 12 the boundary between the regions II and III is given by the 
z = ±oo line. This change of the sign of the variables z (and t) is not troublesome for the continuous extrapolation 
of Eq. (Id. 1511 . because we obtain 


V'out(z,w)exp(ia;log|t|) ~ '0in(z, w) exp(*a; log |t|) ~ exp(iwr), 


(3.16) 


in the limit z ^ ±cx) where we use V —>■ ±cx) in Eq. (ITTIl . (The amplitudes gont{z, lu) and giniz, to) are also continuous 
at the boundary z = ±cxd.) Hence, Eq. 113.1511 holds in the region III without any modifications, except that the 
integral in the second term of the right-hand side should be understood as the contribution from the source point 
{s,y) located in the regions II and III. 

To clarify the physical properties of the solution a{t,z), it will be convenient to use the null coordinates u and v 
for the observation point and u' and v' for the source point. This coordinate transformation gives the Jacobian 


d{s,y) 
d{u', v') 


J{u’, v') 


sy{l-V^{y)) 

2V{y)u'v' 


(3.17) 


and Eqs. iimii and 113.1511 are rewritten into the unified form 


a{u, v) 



dv'H{u, v\u\ v')Qi{u', v')J{u', v'), 


(3.18) 


where the lower limit of the integral with respect to v' should be u', because the regular center y = 0 corresponds to 
v' = u'. If the observation point {u,v) is located in the region I, the function H{u,v\u',v') in Eq. 113.1811 is equal to 
the original Green’s function G(t, z\s,y). However, if the observation point (u, v) is located in the region H (or in the 
region HI), we obtain 


H{u,v\u ,v') = Gi{t,z\s,y), (3.19) 

for the source point {u',v') located in the region I, and 

H{u,v\u,v') = Gi{t,z\s,y) - G 2 {t,z\s,y), (3.20) 

for the source point (n', v') located in the region H (or in the regions H and HI). The validity of the integral represen¬ 
tation 18|1 is supported in Appendix El by applying this formula to a static held a{t, z) = a{r) in the Minkowski 
background. 

Our main concern in this paper is the emission of outgoing radiation toward the future null inhnity. Hence, using 
the two modes ipo^tizjto) and ^in(z,a;) instead of '0out(z:, w) and '0in(z,a;), we divide the frequency-domain Green’s 
function given by Eq. (E3 into the outgoing and ingoing wave parts as follows. 


Gj{z,y,uj) = tfjont{z,uj)G^out{y,oj) + iloin{z,o;)Gjin{y,uj), (3.21) 

where 7 = 1,2. The outgoing part of the right-hand side of Eq. (liOTl) can derive the “retarded” Green’s functions 
Gjoutiu,v\u',v') as 


G^qh^(^u^v\u ,n ) 


1 

27rs 





which reduces to 


G.yOUt(^|^ ) 


1 

27rs 


gout(^2,a;)G^out(2/,a;)e*“'°8(-“/l*l)da;, 


(3.22) 


(3.23) 


in the limit v —>■ 00 (i.e., z ^ Z 2 ) with a hxed u. Via the above-mentioned procedure and Eq. we can 

obtain the outgoing part aout(u,v) of the vector potential a(u,v). The hnal formula for the outgoing radiation held 
aout(w) = aout(M, Z 2 ) observed at the future null inhnity v ^ 00 (i.e., z ^ Z 2 ) is given by 


O-out {tt) 




dv'Houtiu\u', v')Qi{u’, v')J{u\ v'), 


(3.24) 






where Ho^t{u\u',v') is equal to Giout{u\u',v') for the source point {u',v') located in the region I, and is equal to 
Giout{u\u', v') — G 2 ont{u\u'^ v') for the source point («', v') located in the regions II and III. The explicit form may be 
written by 




I 




gout (-2:2, 

w{uj) 




,iui log(-«/|u'|) 


■ N{y,uj)e^ 


i\og(u/- 




duj. 


(3.25) 


where 


K = 0 , N = {bi-b2)9in{y,co), (3.26) 

in the range v' > —u' (i.e., in the region III). If the source point is located in the region I or II, the functions K and 
N should be expressed by (7out(g,w) and gi-a^{y,uj) instead of gout(gj<x;) and gi^{y,uj) to be regular at the boundary 
y = zi (i.e., u' = 0). Hence, we obtain 


K = {bi- I)gout(g, w), 

in the range 0 < u' < —u' (i.e., in the region II), and 


N = bigin{y,uj), 


K = 6 igout(g,w), 


N = 6 igi„(g,w), 


(3.27) 

(3.28) 


in the range u' < v' < 0 (i.e., in the region I). It should be also remarked that from the causality condition we can 
use Eq. only in the range u' < u for the retarded time u' of the source point. 

Because the physical quantity to be observable is the energy flux of the outgoing radiation, we must calculate the 
derivative of a{u) with respect to the retarded time u. Here, for later convenience, we denote the derivative with 
respect to log(—u) by a dot, and we have 

dout(M) ='I'(u) + <&(«), (3.29) 

where 

/ U pOO 

du' / dv'HoMt(u\u',v')Qi{u',v')J{u',v'), (3.30) 

-00 J u' 


and 

poo 

$(u) = u / dv'Hont{u\u,v')Qi{u,v')J{u,v'). (3.31) 

J U 

The future Cauchy horizon due to naked singularity formation appears on the boundary z = Z 2 (i-e., at the retarded 
time u = 0) of the region HI. Hence, to understand clearly a distinctive effect of the naked singularity, our task in the 
following sections is to investigate the asymptotic behavior of Oout (w) in the limit u —> 0 for clarifying whether or not 
a burst-type emission of the infinite energy flux can occur in the self-similar collapse. 


IV. DECOMPOSITION OF CONTRIBUTIONS TO THE GREEN’S FUNCTION 

Now let us try to obtain explicitly the asymptotic behavior of the functions 'I'(u) and d)(M) in the limit u ^ 0 
according to the formalism given in the previous section. The first step would be to calculate the integral with respect 
to w in Eq. (E 2 S 1 ), by considering the analytic structure of the Wronskian factor w in Eq. as a function of w. 
Hence, following the analysis of test fields in black hole background, we decompose the Green’s function into several 
distinct parts corresponding to contributions from high-frequency waves with |a;| —> 00 and “quasi-normal” modes of 
the system with complex frequencies giving w = 0. 


A. High-frequency contribution 

In the high-frequency limit |a;| ^ 00 we can apply the WKB approximation to obtain the homogeneous solutions 
for Eq. (13.411 . and at the leading order the amplitudes gout(z,<^) and gin(z,uj) are found to become constants denoted 
here by Cout and Cin, respectively, without a mode mixing (i.e., bi ~ I, 62 ~ 0). Then, the Wronskian factor w is 
approximately given by 


w{uj) ~ 2ia;CoutCin 


( 4 . 1 ) 
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FIG. 3: The source region shown in the Penrose diagram, which is given by the null coordinates u and v'. The null geodesics 
v' = u and u' = u starting from the past null infinity and arriving at the observation point denoted by (u, Z 2 ) are shown by 
the thick solid lines. The broken line corresponds to the v' = —u line. The boundaries of the source region A are given by 
u' = — 00 , v' = u, u' = u and v' = —u. 


It becomes easy to calculate the derivative of Hout with respect to log(—u) (denoted by Hout) under this approxima¬ 
tion, and we have 




1 


25g!v(y)-A(y) I g. 




(4.2) 


which is applied to the estimation of II'(m) in Eq. 11.3.3011 . Denoting the high-frequency part of by 'Ff, we arrive at 
the result such that 'I'f(u) = T'fi(u) + '^F 2 (u), where 


'I'fi(u) = 


Cout^ 

2 Cin . 


Qi(u',u)J(u',u)du', 41^2 (m) = 


Qi(u,v')J(u,v')dv', 


(4.3) 


with the integrand factor Qi defined by Qi = Qijse‘'~^ = for a given current distribution ji. 

High-frequency waves will be able to propagate without scattering due to the spacetime curvature and the centrifugal 
barrier (depending on Z). In fact, the i5-functions which appear in Eq. 14.211 allow us to use the approximation of 
geometrical optics for the propagation of high-frequency waves to the future null infinity at the retarded time u (along 
the path drawn in Eig. 0). The source points {u', v') giving the first term d^Fi in 41^ are located on the ingoing null 
line v' = u. Hence, the outgoing flux represented by 4 >fi is a result of the reflection (at the regular center v' = u' = u) 
of ingoing null rays propagating along v' = u, which are generated by the current distribution ji at the source points 
{u',u) ranging from u' —00 to u' = u. On the other hand the second term ^^2 represents the outgoing flux which 
is generated by the current distribution ji just on the past null cone u' = u (ranging from v' = u to v' ^ 00 ) and 
is directly propagating to the observation point. If the retarded time u of observation at the future null infinity is 
nearly equal to zero, the source points for 41^1 (or 41 ^ 2 ) on the ingoing null line v' = u (or the outgoing one u' = u) 
become nearly identical with the past (or future) Cauchy horizon, where we obtain roughly m' ~ s (or v' ~ s) except 
an unimportant constant factor. Then, the Jacobians J in the integrands of Eq. are roughly given by 


J(u,u)~-(-^) , J{u,v') , 

u \w / u \ v' / 

where ki and K 2 are the values previously defined in Sec. m Hence we obtain 

pu pOO 

4'fi(m)(-■ u)'^W ji{u',zi){-u'y~‘'^^du', 4'f2(m) ^ (-m)'‘M ji{v', Z 2 )v''^~'"^dv'. 


(4.4) 


(4.5) 


^ —00 J U 

As was mentioned in Sec. QH we assume a rapid decrease of the current distribution ji near the regular center. A 
plausible choice would be that ji(s,y) ~ h(s)r‘^{s,y)S‘^(]j) at the source point r = sy ~ 0, where h{s) is introduced 
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as some amplification of ji in the process of the self-similar collapse. It may be interesting to consider the case that 
ji{s,y) diverges at the onset s = 0 of the naked singularity for a fixed nonzero y (namely, hs^ —> oo as s —> 0). 
However, our purpose in this paper is to pursue the possibility of amplified generation of electromagnetic radiation 
due to the self-similar dynamics of background geometry involving a naked singularity, which requires us to calculate 
'ki?i and f 2 without any divergence of the source current jf. (Hereafter we assume h to be finite at s = 0 for 
simplifying our discussion.) Furthermore, we obtain the relations 

Kl=l -^-, K2 = 1 -(4.6) 

where Ryv and Ruu are the components of the Ricci tensor evaluated at z = zi and z = z^, respectively. Then, if the 
strong energy condition holds for collapsing matter, we can require the inequality 

1^1 <1, < 1, (4.7) 

to show that the integrals in Eq. 14.511 converge even in the limit u 0. (Of course, it should be also assumed that the 
current distribution j) rapidly decreases as the source point becomes close to the past or future null infinity u' —>■ —oo 
or v' —>■ oo). Hence, finally we find the asymptotic power-law behaviors of 'kpi and ']^f 2 giving the outgoing flux of 
high-frequency waves to be 


4'fi(m) ~ (-«)''% ^'f2(w) ~ (-w)'^^ 


(4.8) 


in the limit u 0. 

Here let us turn our attention to the additional term <!)(«) in Eq. 15.2911 . Note that the source points (u',v') giving 
<!)(«) are also located on the past null cone u' = u in the same way as 'I'i^’2 in Eq. 14.511 for high-frequency waves. 
The contribution to $(u) given at the retarded time u will represent outgoing waves generated on the past null cone 
ranging from v' = u to v' oo. This allows us to estimate easily the asymptotic behavior of $(u) in the limit m ^ 0 
without using the high-frequency approximation. The key point is that in the limit u 0 the range of the past null 
cone in the regions I and H shrinks to zero (see Fig. 0 ). Hence, using Eq. valid in the region HI, the function 

ffout('ulu,v') in Eq. 15.5111 can be given by 


ffout(ulu,v') 


I 

27rse‘"(^2)-A(22) ’ 



(&i - b 2 ) 
w{uj) 


gout{z2,0j)gin{z2,uj)duj. 


(4.9) 


in the limit u —*■ 0, where s v' and y = Z 2 for the source point u' = 0. The convergence of the integral I with 
respect to uj in Eq. 63 may be subtle, because the integrand is proportional to 1 /uj in the high-frequency limit. 
However, the contribution of such a high-frequency part can be canceled out, for example, if the infinite integral of 
any function cr(w) is defined as 


/ CX) pUJQ 

a{u!)duj = lim / a(uj)du 

-oo 


By virtue of such a regularization of the integral we obtain 

21 


$(-«) = -4'i7’2(w) ~ (-^t)'"^ 

TT 


(4.10) 


(4.11) 


in the limit u —> 0. If wave frequencies are not so large, outgoing waves will be efficiently back-scattered, and their 
contribution to the function Hautiu\u,v') in <I>(m) becomes important to estimate the value of I. Nevertheless, from 
Eq. (ICTTl we can conclude that such a back-scattering effect is not relevant to the asymptotic power-law behavior of 
$(m) which becomes identical with that of '^f 2 {u). 


B. Contribution from poles in the Green’s function 

The frequency-domain Green’s function may contain singularities in the complex w-plane. In black hole background 
it is well-known that late-time tails of massless test fields are generated by a low-frequency contribution to the 
frequency-domain Green’s function, which appears as a w = 0 branch cut in the Wronskian factor w{uj). There may 
be also a branch cut in the frequency-domain Green’s function 63 considered here. However, in this paper we do not 
pursue such a problem. Here our investigation is focused on the pole contribution to '^(u), which will be interesting as 
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a resonant behavior of test fields in the self-similar background. We use the residue calculus to evaluate the integral 
(with respect to ui) in i?out, from which the high-frequency part is subtracted. 

Recall that the simple poles in the frequency-domain Green’s function are given by the zeros of w{uj). We must 
remark that even in the Minkowski background there are zeros at w = ±* (see Appendix and the two mode 
functions ^out(z,±*) and ipin{z,±i) given at the frequencies can satisfy the relation tpQ = i/^out + V'in = 0. (There is 
also a zero of ic at w = 0, but it is not a pole because the factor w in w is canceled out by virtue of the differentiation 
of Eq. Id.2511 with respect to log(—rt).) As shown in Appendix IBI through the analysis of '0out and i/jjn near z = Zi, 
if the background spacetime is dynamical, the zeros of w{uj) should change to u) = ±Kii (note that the value of ki 
becomes unity in the limit to the Minkowski background). From the same analysis of the different pair given by 
Eq. (HTT^ near 2 = ^2 h is easy to see that the outgoing mode '(/'out remains linearly independent of the ingoing mode 
ipiri at the frequencies to = ±kij, because in general the equality ki = K 2 does not hold. Hence, the requirement 
such that 0 = 'ipo = ipo^t + '0in = V'out + V'in = 0 at w = zLkii can become consistent only if both gout(^, ±Kii) and 
gin{z,±Kii) vanish at any z (in the same way as ipo)- It is also required that the coefficients bi and 62 in Eq. (Id.fill 
must be inversely proportional to the factor -I- for keeping the functions goutiz^uj) and gin(z,<x>) finite even at 
u! = ±Kii (namely, foi^out and 625111 become finite in the limit uj ±Kif). Then, the residue calculus of Eq. (Id. 2 511 
for the zeros of w at w = ±Kii leads to the result that the pole contributions (denoted by Hott and i^out) lo ^out 
exist only for the source point located in the range |r’'| < |u| < \u'\ (namely, in the region A drawn in Fig. 0 ), and 
they are written by 




Kii lim^ {61 {io) ffout {z2,(zi)} 5out (5, ±Kii) 


—u 




(4.12) 


where w' = dw/duj, lim^ means the limit w —*■ ±kiz, and the relation 


5out iy,±Kii) = -gin{y,±Kii) 


(4.13) 


is used. 

Though the physical interpretation of the modes satisfying the relation tpout = ~'0in at w = ±Kii is unclear, their 
contribution giving represents a prompt emission of outgoing waves in the region A including the past Cauchy 
horizon. In particular, for the growing mode with Im(w) = > 0 the function (cx can infinitely increase 

as the retarded time u of the observation point approaches zero. However, the source region in the range u < v' < —u 
is restricted to the past Cauchy horizon in the limit u —> 0, and the cu = kA contribution (denoted by 4'+) to 4 '(m) 
in Eq. (TOnil turns out to have the asymptotic behavior 

4/+(n) -(_„)-!, (4.14) 


because the Jacobian is given by J{u',v') ~ (l/n')(n'/u')”i near v' = 0, and 5out(y, «!*) is finite at 5 = zi (see 
Appendix 01). This is quite similar to the high-frequency contribution which is a result of the reflection of an 
ingoing flux (generated on the null line v' = u) just at the regular center u' = v' = u. The u = kA contribution also 
may be a result of the conversion of an ingoing flux (generated in the region A) into an outgoing flux, which occurs 
on the null line u' = u with the finite range u < v' < —u. 

In the calculation of the contribution (denoted by Jt-) from the decaying mode with Im(a;) = —ki < 0, we must 
remark that 5in(y, ~s^i*) is finite aX y = Zi (see Appendix 01, thus from Eq. 114.131) 5out(5 ,-kr) ~ \v' near 
v' = 0. Then, we find the asymptotic behavior such that 4'_(u) ^ (—in the limit u —> 0. It is obvious that this 
oj = —kA contribution becomes unimportant in comparison with 

In the above-mentioned analysis we have revealed the contributions from various source points located in the range 
M < n' < 00 to the outgoing flux observed at the future null infinity. It is shown here that we can also obtain a 
contribution from the source point located in the range u' < v' < u corresponding to the inner part of the region I, 
if there exist quasi-normal modes of the dynamical self-similar system with complex frequencies giving w(ui) = 0, at 
which the function ipout{z:,uj) outgoing at z = zi becomes regular at the regular center z = 0. Namely, we assume the 
relation Ain = Ao — Aont = 0 aX to = ujq (5 = 1, 2, • • •), instead of the relation Ao = 0 at to = XzkA- 

The real part of ujq for such quasi-normal modes may be nonzero, and it is proved in Appendix0]that the imaginary 
part of ujq should be negative. Further, from Eqs. (IXTTl) and (imi) . the condition Ain = 0 leads to the result such 
that 5i = 62 = 1- Then, using Eq. and the residue calculus associated with these poles, we obtain the uj = ujq 

contribution (denoted by iJout) fo iJout as follows. 


Hont(.u\u\v') 


^q9out{z2-, ^q)9outiy•! ^q) 



X 


(4.15) 
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Note that Eq. (14.1511 is valid only in the range u' < v' < u, and Hout = 0 at any other source points in the regions 
I, II and III. The quasi-normal modes represent a disturbance excited in the inner region around the regular center, 
which can generate outgoing waves arriving at the future null infinity. 

The key point to calculate the contribution (denoted by 4',) from the quasi-normal modes to 'I'(u) is that the 
integral in Eq. iionii with respect to v' is limited to the range u' < v' < u, and the integrand becomes proportional 
to the factor in the limit v' 0. Then, it is easy to see that the asymptotic behavior of in the limit 

M ^ 0 is given by 

(4.16) 

if the quasi-normal mode decays slowly (namely, |Im(a;g)| < ki). Interestingly, we find an oscillatory behavior 
'kq ^ iog(-“) fQj- the outgoing flux observed at the future null infinity. However, if it decays rapidly (namely, 

|Im(a;g)| > ki), we obtain again the power-law behavior given by 

4',(u) ~ (-u)"E (4.17) 

Then, no new feature appears for the time evolution of 'I'(u), even if the quasi-normal mode is generated. 

In summary, we have obtained various contributions to dout(w) representing the outgoing flux observed at the 
retarded time u in terms of the decomposition of the Green’s function as follows, 

dout = «'f 2 + «> + ^Fi + ^'++ XI (4-18) 

where the unimportant contribution is neglected. According to the difference of the source region, they show 
different asymptotic behaviors in the limit u 0. The first two terms f 2 and <I> due to outgoing waves originated 
on the past null cone u' = u show the power-law behavior with the power index K 2 , while the next two terms 4'^?! and 
with the source region near the past Cauchy horizon v' = 0 show the power-law behavior with the power index 
Ki. If slowly decaying quasi-normal modes are excited near the regular center, an oscillatory behavior with respect to 
the logarithmic time log |m| can appear in the final term 'kg. The important point is whether or not these terms 
can induce a divergent energy flux in the limit u ^ 0, at which the effect of formation of a naked singularity appears. 
In the next section we will define the relation of dout to the outgoing energy flux measured by distant observers and 
discuss the values of K 2 , ki and ujq under specific self-similar models. 


V. THE OUTGOING ENERGY FLUX 


In the previous sections we have used the comoving null coordinates u and v for specifying the observation point. 
The metric component g^v for the double null coordinates is written by 


9uv — 


uv 


1 ) 


(5.1) 


The comoving coordinates become crucially different from the double null coordinates u and v used by a “static” 
observer present at a point sufficiently distant from the center, where we obtain ~ 1. As was shown in 0 , if the 
spacetime describing a self-similar collapse is smoothly matched to the Schwarzschild spacetime at the star’s surface, 
the retarded (or advanced) time u (or v) becomes equal to the retarded (or advanced) Eddington-Finkelstein null 
coordinate near it = 0 (or near i;=0). Hence, in this section we evaluate the outgoing energy flux F{u) per a unit 
time of It, instead of it, which is important as a quantity measured by the static observer near the future null infinity. 
Near the future null infinity (namely, at z ~ Z2), Eq. leads to guv ^ {\u\/vY'^ and we find the relation 


|ii| ~ lurT (5.2) 

Hence the outgoing energy flux F(u) is given by 

/ dUout \ f dUout \ f- 

~(fl) 

It is also remarkable that near the past null infinity corresponding to z ~ Zi, we obtain |i;| ~ |ii|"’'. 

It is clear from Eq. that the outgoing flux defined by doout/dit diverges in the limit u ^ 0, because in general 

the parameters ki and K 2 become smaller than unity, as was mentioned in Eq. (ITO . For example, the first two terms 
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'I'F2 and d) in Eq. iBTEl give the divergent term |u|”^ ^ to duout/du, which will represent a result of prompt emission 
of outgoing waves (measured by the comoving time u) on the past null cone u' = u in the range v' > u, where the 
coordinates u' and v' denote a source point for wave generation. However, the ratio u/u of the two retarded times in 
Eq. 15.311 works as a red-shift factor in calculation of the energy flux defined by F{u), and the d>F2 and $ contribution 
(denoted by fi{u)) to duout/du remains finite even in the limit u ^ 0 . 

On the other hand the next two terms 'I'fi and in Eq. (OSI) give the divergent term ^ to daout/du. As was 

mentioned in the previous section, these terms will represent a result of the reflection of an ingoing flux propagating 
along the null line v' = u ^ 0, where we obtain |i;'| ~ for the advanced time v' of the static observer. Thus, the 
divergent term may be interpreted as a blue-shift factor v' jv' for ingoing waves propagating along the path 

v' = u toward the center. Of course, the red-shift effect also becomes important when the reflected waves propagate 
to the future null infinity, and the 4 'fi and contribution (denoted by /2) to duo^t/du is given by 


h{u) 


u 

u ^ |u| 





(5.4) 


in the limit u 0. If K 2 is smaller than ki, the contribution /2 becomes zero rather than finite as ft —> 0, which 
contrasts with the limit of fi. This may mean that the red-shift of the reflected waves is more pronounced than that 
of the waves without passing through the center. 

Let us denote wave frequencies for quasi-normal modes by ujq = ojRq iujiq for real ujRq and ujiq. If \ujiq\ is larger 
than Ki for all q, the 'kg contribution (denoted by fsq) to duout/du shows the same power-law behavior as /2. We 
note that the amplification due to the blue-shift effect on the null line v' = u can still work even for rapidly decaying 
quasi-normal modes. If \u}iq\ becomes smaller than ki for some q, the contribution f^q in the limit u 0 can be 
written by the oscillatory form 


/3,(u) ~ lul^’cos ("^^loglulV ^ =1^^-I. (5.5) 

\ K2 J K2 

Owing to prompt generation of waves induced by slowly decaying quasi-normal modes in the range u' < v' < u, the 
blue-shift effect on the null line v' = u becomes insignificant. 

In conclusion, we find the asymptotic behavior (in the limit u —> 0) of the outgoing energy flux F{u) measured by 
static observers present near the future null infinity as follows, 

F{u) (X ^fiiu) + f 2 iu)+ '^f 3 q{u)^ - |ci-l-C 2 |u|^°-|-^C3|?2|^‘>cos^^log|M|^| , (5.6) 

where ci, C2 and C3 are constants. The third term in Eq. (1^ appears only if there exist quasi-normal modes satisfying 
the condition \uJiq\ < ki. By virtue of the red-shift effect no divergence of the energy flux F{u) occurs even at the 
moment when the effect of naked singularity formation arrives at the future null infinity, unless ki or \u}iq\ becomes 
smaller than K 2 - For the self-similar collapse such that ki becomes smaller than \u}iq\ and K 2 , the power-law divergence 


F(m) - (5.7) 

will be observed as a dominant behavior with the power index limited to the range —2 < 2^o < 0- If slowly decaying 
quasi-normal modes are excited (namely, if is smaller than ki and K2), tbe oscillatory divergence 

F{u) ~ cos^ jQg ^ _2 ^ < 0, (5.8) 

can become a more interesting precursor of naked singularity formation. Of course, the generation of a divergent 
energy flux means that the naked singularity formation is an unstable process against a backreaction of perturbations. 
Nevertheless the above-mentioned precursor phenomenon may significantly appear, before the backreaction effect 
becomes important. Even if F(u) remains finite at it = 0, the anomalous time evolution of F(u) due to the damping 
term or cos(a;_R5 log |m|//C 2) will appear in principle as an observable effect different from usual black hole 
formation. 

Here we briefly comment on the outgoing energy flux in black hole formation in terms of Eq. lEH). This corresponds 
to the case V'(z 2 ) 0 (namely, Zm = Z 2 = Z 3 in Fig. for which the Cauchy horizon coincides with the event 

horizon. Then, the power indices ^0 and ^q become infinitely large, and the power law decay of the second and third 
terms in Eq. 15.61) will change to an exponential decay. The first term ci in Eq. 15.611 also vanishes in the limit 
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^0 



FIG. 4: Numerical evaluation of the power law index for dust collapse parameterized by k. The value of diverges at 
k = k*, where the future Cauchy horizon coincides with the event horizon. 


^'{^ 2 ) 0, because the Jacobian J{u,v') involved in the integrals (I4.dll and Id.dill (giving 'I>i7’2 and 4)) is found to 

be proportional to the parameter V'(z 2 ) at the leading order in the limit u —> 0, by considering the constant factor 
omitted in Eq. 14.411 . The formula given by Eq. (1^ can describe the infinite red-shift effect at the event horizon 
which makes the outgoing energy flux F(u) completely vanish at it = 0. 

The key parameters for determining the asymptotic evolution of F(u) are the ratios ki/k 2 and \ujiq\/K 2 - To discuss 
the values of the parameters, we must specify the self-similar model. If the collapsing matter is a perfect fluid, the 
self-similarity (1^ is satisfied only for the equation of state P = ap, and we find the inequality ki > K 2 for 0 < a < 1. 
The general proof of this inequality except the case a = 0 is given in Appendix |0 For the collapse of dust (namely, 
for a = 0) we can obtain explicitly the metric 12.111 as follows. 



from which the function V(z) is derived as 

If Qz 1 

where k is a positive constant. The global structure described in Fig. El (corresponding to the behavior of V shown 
in Fig. ^1 is possible only in the range k > k* = (26 -I- 15-\/3)/6 ~ 8.7, and from Fig. 0it is easy to see that the 
power law index ^0 becomes always positive. Naked singularity formation due to self-similar collapse of a perfect fluid 
will not induce a divergent outgoing energy flux F(m), unless there exist slowly decaying quasi-normal modes with 
< 1- It will be interesting to check whether or not the inequality ki > K 2 holds for any other collapsing 
matter (for example, in the case of a scalar-field collapse). 

To analyze the quasi-normal modes in general, it will be useful to rewrite the homogeneous equation for the outgoing 
mode ^out(z,w) into an one-dimensional Schrodinger-type differential equation with a potential U{z), which is given 
by the metric components iz, S and A and decreases in proportion to 1 — V^ in the limit z —>• zi (see AnnendixIBlfor the 
equation). The existence of the quasi-normal solution ipout satisfying the regularity at z = 0 is shown in Appendix I dI 
by assuming some suitable form of U. We find the case such that the absolute value \uJiq\ of the imaginary part of the 
complex frequency coq with nonzero LOuq can be smaller than ni. However, it remains unclear whether or not such a 
case is possible for the metric satisfying the Einstein equations. This problem would be more extensively investigated 
in future works. 

We have developed the Green’s function technique to study time evolution of electromagnetic perturbations in 
the self-similar background, which will be applicable to any other test fields (for example, scalar and gravitational 
fields). We have also clarified that the various contributions to the Green’s function (for example, high frequency 
part and quasi-normal modes) given by Eq. (14.1811 are closely related to the difference of the source region where the 
corresponding perturbations are generated. 
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Our main result is the finding of the key parameters k,iIk 2 and \uJiq\/K 2 to determine the occurrence of divergent 
energy flux associated with naked singularity formation. The value of the complex quasi-normal frequency Wq should 
depend not only on the background geometry but on the test field (and I parameterizing the field angular momentum). 
On the other hand, the values of ki and K 2 are purely geometrical quantities given by the derivative of the metric 
on the past and future Cauchy horizons, respectively. The blue-shift and red-shift effects defined hy v/v and u/u for 
ingoing and outgoing waves propagating near the Cauchy horizons are parameterized by these geometrical quantities. 
The competition of such two effects in generation of outgoing energy flux would be an important feature of naked 
singularity. 

It is interesting to compare our result with the result giving by Chandrasekhar and Hartlel^^, who have considered 
time evolution of electromagnetic and gravitational perturbations inside the Reissner-Nordstrom black hole. The black 
hole contains a timelike singularity visible by a falling observer in the black hole. Estimating the energy flux received 
by an observer near the (future) Cauchy horizon for the timelike singularity, they have found that its amplitude has an 
exponential behavior with the rate k_ — k_|_, where k_|_ and k_ are the surface gravity factors evaluated on the event 
horizon and the Cauchy horizon, respectively (note that the event horizon corresponds to the past Cauchy horizon). 
Because the rate k- — k+ is always positive, the energy flux diverges when the observer arrives at the Cauchy horizon. 
The power index which is written by the difference of the constant k evaluated on the future and past Cauchy 
horizon by Eq. (EUl), resembles the rate, except that it is not always positive. This commonality would provide an 
important insight into an universal feature of perturbation response in various geometries with a naked singularity. 

Eurther we would like to emphasize that the power index 7 of the divergent energy flux F ^ with (or without) 
an oscillation (corresponding to 7 = 2^q (or 7 = 2^o)) is allowed in the range — 2 < 7 < 0. This sharply contrasts 
with the semiclassical result giving the unique power index 7 = —2 for the radiation of quantized test fields without 
depending on the parameters of the self-similar background [ill lia ll^ Il4 IlM . Such a difference of the power-law 
divergence may become an important clue for understanding profoundly quantum properties of perturbations around 
naked singularity. 


APPENDIX A: APPLICATION TO MINKOWSKI BACKGROUND 


Let us apply Eq. (ixra to Minkowski background iz = A = 0 and S' = 1 for checking the validity of the formulas 
given in Sec. imi Here we consider the dipole field with 1 = 1. Then, the homogeneous solutions for Eq. (lO) are 
easily given by 


% 2 

gout(2,w)=w-, gin{z,oj) =uj + (Al) 

Z 2: 

Of course we obtain 61 = 62 = 1 for the coefficients in Eq. lITTTll . and the Wronskian factor is given by 

w{ijj) = 2iu!{uj + i){uj — i). (A2) 

Then, from Eqs. (in . lIXT^ and we find 


H{u, v\u', v') 


2{uv + u'v') — {v' + u'){v -I- u) 
{v — u){v'— u'){v'+ u') ’ 


(A3) 


for u < v', and H{u,v\u\v’) = 0 for v' < u. To calculate explicitly Eq. 13.18|l . we assume the shell-like current 
distribution ji{t,z) = joS{r — rg) for a constant jo- Then, we find a(r) = —47rjoro/3r for r > rg, and a(r) = 
—47rjor^/3ro for r <rg. This corresponds to the vector potential for a static dipole magnetic field. 


APPENDIX B: ANALYSIS OF HOMOGENEOUS SOLUTIONS 

Here we analyze the homogeneous solutions 'ip{z,uj) for Eq. (111.411 . by rewriting it into the one-dimensional 
Schrodinger-type differential equation 


^ - Uix)) ^ = 0, 


where the new variable x and the “potential” U are given by 

1 


^ = 2 i^ont{z) - hiniz)) , 


(Bl) 


(B2) 



and 


16 


U{x) = 




Note that t7 > 0 in the region I shown in Fig. [3 The function ip is defined by 

iP{x,u) = 


(B3) 


(B4) 


To obtain the homogeneous solutions giving ipo and ipout, we must consider the behavior at the points z = 0 and 
z = zi corresponding to cc = 0 and x = oo, respectively. The local flatness at the regular center and the suitable gauge 
choice allow us, without loss of generality, to assume that iz ~ 0 and ~ mS ~ C'a(— near z = 0^, where C\ 
and m are positive constants. Then, we find the approximate form 


U{x) 


l{l + l) 


near a; = 0. On the other hand it is easy to see that 


(B5) 


U{x) ~ (B6) 

in the limit x oo, where Uoo is a positive constant. It is remarkable that we have U = l{l + l)/sinh^a; and 
z = — tanha; for the Minkowski metric. 

Denoting the solutions for Eq. (IB 111 by '0out and tpin corresponding to '0out and '0in, respectively, we obtain 


ipont{x,uj) = gout{z,uj)e'^‘^^, '0in(a;,u;) = gin(z,a;)e (B7) 

The approximate forms of these functions for large x are given by 

5out(a;,a;) ~ g^^iuj) |l + . (B8) 

and 

gUx,u;) o. g-M |l + , (B9) 

for some functions g^f. and dependent on uj. To keep these functions well-defined even at a; = ±Kii, the functions 
and <7;^ are required to contain the factor Ki — iuj and Ki+iuj, respectively. Then, we obtain w(w) = 2iuig'^^g'^ = 0 
at a; = ±Kii, for which it is easy to see that the two modes ■i/'out and V’in become linearly dependent at the leading 
order in the limit x —>■ oo. 

Next, using Eq. m. we prove that the imaginary part of the frequencies ojq of quasi-normal modes denoted by 
ipq{x) should be negative. We require the boundary conditions for the quasi-normal modes such that ipq ~ x^'^^ near 
X = 0, and ipq ~ in the limit x —*■ oo. Then, integrating Eq. m multiplied by the complex conjugate ip* in the 
range 0 < x < xq, where xq is a positive constant satisfying ipq{xo,ujq) ^ 0, for sufficiently large xg we obtain 


Xuj^ + iujqV — Z = 0, 


(BIO) 


where the coefficients given by 


X = 




dx, 


Y = 


^q{xo) 


z = 


dipq 


dx 


U 


Ipq 


dx, 


(Bll) 

are positive definite. It is easy to see that Im(wq) = —Yj2X < 0 for AXZ > and Im(a;g) = —{Y ± 
VE2 - 4:XZ)/2X < 0 for AXZ < 


^ Then, via the coordinate transformation r = we can find that near z = 0, 


ds^ ~ —dt^ + 




(df"^ +f'^dn^). 


m 
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FIG. 5: Schematic description (sold lines) of the functions V+(z) and V-{z). The dotted line denotes the function V{z) for 
z < 0. The line V = V-{z) is the point symmetric line to the dotted line with respect to the origin. 


APPENDIX C: PERFECT FLUID COLLAPSE WITH PRESSURE 

In this appendix, we prove that the value of K 2 is always smaller than that of ki in the spacetime describing the 
collapse of a perfect fluid with the equation of state P = ap for 0 < a < 1. For convenience, we express the function 
V{z) for z > 0 and the point symmetric function to the function V (z) for z < 0 concerning the origin (z, V) = (0, 0) 
by the notation V+(z) and V-{z), respectively (see Fig. EJ. As stated in the function V-{z) is also given by the 
metric satisfying the Einstein equations. In addition, we define' = zd/dz only in this appendix. Hence, what we show 
is that the value of V+(z2) is always smaller than that of U_(|zi|). 

Now we start the proof with a similar notation used by Carr and Coley in |^. They introduced the functions A(z) 
and B{z) defined by (47rpr^)““/(^+“^ = and S{z) = for positive constants Aq and Bq. From the 

Einstein equations, the functions I'iz) and A(z) can be expressed in terms of the functions A(z) and B{z). Hence, the 
function V(z) for z > 0 can be expressed as 

U(Z) = 

for a positive constant Vq (see Eq. (2.19) of [13). In addition, we use one of the 

which is shown in Eq. (4.23) of [13. When V = 1, Eqs. IICl ll and IIC2ll give 

V =1- 

1 -|- ct 

This means that the values of V-{z) and V^{z) at 2 = |2i| and 2 = 22, respectively, become 

I/_ (|zi|) = 1 - ^^|zi|-(i-“)/(i+“)5-2(zi), (C4) 

\/o(l + a) 

and 

^+(^2) = 1 - —(C5) 
Vb(l + a) 

As stated in [13, the solution describes the monotonical collapse. This means that the value of S'(zi) is larger than 
that of S'(z2). Hence, when a = 1, we find the inequality U+(z2) < V- (|zi|) from Eqs. (EH) and ESI). On the other 


(Cl) 


Einstein equations 


- 1 




(C2) 
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hand, for 0 < a < 1, Eqs. inn) and insi) indicate that if the value \zi\ is larger than the value Z 2 , the value of 
V^{z 2 ) is always smaller than that of V-{\zi\). The condition \zi\ > Z 2 is satisfied if the inequality V^(z) > V-{z) 
holds for the whole region z > |zi|, as depicted in Fig. [3 Therefore, the remaining problem is to show the inequality 
V+(z) > V-{z) in the whole region z > |zi| for 0 < a < 1. 

We firstly show that V+(z) > V-{z) holds at the sufficiently large z. By requiring the left hand side of Eq. (EH) to 
be finite as E ^ oo, we find that at the sufficiently large z, Eq. becomes 


B = 


A 


which is also shown in Eq. (4.25) of [^. From this equation and Eq. m, we can obtain 


(C6) 


V 

V 


1 — a 


1 + a 


- 2aB. 


(C7) 


Hence, we find 


VI v; 

V-~V^ 


-2a(B'_ - B:^), 


(C8) 


where the functions B-{z) and B+{z) denote the values of B{z) along the lines of V{z) = E_(z) and V{z) = V+{z), 
respectively. The monotonical collapse also means that B'_ < 0 and B’^ > 0. In addition, the equality V+{z) = V_(z) 
holds in the limit z ^ oo. These facts means that the inequality V!_{z) > E(.(z) holds. Hence, the inequality 
E)-(z) > V_(z) holds at the sufficiently large z. 

Next, we show that the lines V = V+{z) and V = E_(z) do not cross each other in the region z > |zi| to complete 
the proof. We eliminate the functions A{z) and A(z) from Eq. EH), using the function V{z) and V{z), and obtain 


- ) ■ (C9) 

1 + a J 

We consider the difference between Eq. (I(19ll expressed in terms of E_ (z) and B_ (z) and Eq. (1(1911 expressed in terms 
of E)_(z) and H_|_(z). Then, for the region E+ > 1 and E_ > 1, we can find the inequality 


1 — a 


(l-V^)B = -^( 
^ ' 1 + a V 




e-2^ - i I - 





> {e! -V^ + (E+e-^^^ - E.e-^^i^j ^ 


(CIO) 


because of the conditions < 0 and > 0. Now we assume that the lines E = E+(z) and E = E_(z) cross each 
other (i.e., E+(z) = E_(z) = Ec > 1) at some points z. At this point, the inequality HCIOII becomes 



^ (^-2B+ „-2B. 


— e 


)• 


(Cll) 


Because of the monotonical collapse condition B-{z) > B+{z), we obtain the inequality V-{z) > E|_(z). However, 
since the inequality E+(z) > E_(z) holds at the sufficiently large z, the derivatives should satisfy E^(z) > E_(z) at 
the point where E+(z) = E_(z). This is a contradiction, which is solved by denying the existence of the point where 
E+(z) = E_(z). Hence, the inequality E-|_(z) > E_(z) holds in the whole region E_(z) > 1. This conclusion ends the 
proof that the inequality E+( 2 : 2 ) < E_(|zi|) always holds for 0 < a < 1. 


APPENDIX D: APPROXIMATE EVALUATION OF FREQUENCIES OF QUASI-NORMAL MODES 

The oscillatory behavior shown by Eq. (ICTi can appear if there exist quasi-normal modes with the frequencies ujq 
satisfying the conditions Re(a;q) ^ 0 and |Im(a;q)| < ki, though it depends on the value of K 2 whether or not the 
divergence of the energy flux occurs. In general it will be a difficult task to find quasi-normal modes as homogeneous 
solutions for Eq. constructed by self-similar metrics satisfying the Einstein equations. Even the existence of such 
modes will be unclear. Hence, in this appendix we consider a simple form of the potential U{x) defined by Eq. (EH). 
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-30 -20 -10 0 10 20 30 


Ke{ujq) Re(a;g) 

FIG. 6: The complex frequencies uiq of quasi-normal modes for the potential U{x) given by Eq. (Dl). The left and right 
panels correspond to the cases Xc = 1 and Xc = 10, respectively, while the other parameters are hxed to be Z = 1, Uoa = 8 and 
Ki = 0.5. 


Our purpose is to illustrate that only a slight modification of U from the potential given by the Minkowski metric is 
sufficient for obtaining quasi-normal modes. 

Taking account of the asymptotic behaviors given by Eqs. IIB5II and IBBII . we can give a plausible approximation of 
Uhy 


U{x) = 


l{l + l) 


C + (Roo-C)(l-e-2«^") 


for X < Xc, 


for X > Xc, 


(Dl) 


4sinh^ (kix) 

where C and Xc are constants. Because the potential should be continuous at x = Xc, the parameter C is given by 

_ / 4,1(1 + 1) sinh^ (kiXc) 


C = e^ 


V 


- C/oo (1 - \ . 


(D2) 


Therefore the potential U given by Eq. m is parameterized by I, Xc, Uoo and ki. 

From the regularity of quasi-normal modes denoted by V’q at a; = 0, we obtain for x < Xc 


■ipqix) = y/xJl+[l/2){ujqX), 


(D3) 


where Ji+(i/ 2 ) denotes the Bessel function. On the other hand, from the asymptotic behavior ijjq ~ obtain 

for X > Xc 


%{x) = gZA^q) (1 - (3,, 132 ,1 - 


(D4) 


where p. = {1 -I- \/4 + (C/ki)}/ 2, and F is the hypergeometric function with j3i and (32 given by 


/3i/32 


U oQ / zto \ 

4kI 


and 


f3i + (32 



(D5) 


(D6) 


By requiring that the functions given by Eqs. 1D3II and fOl) are smoothly matched at a; = a;c> we can find the 
frequencies ujq. The numerical results are illustrated in Fig. We can confirm that there are a number of quasi¬ 
normal modes with Re(a;q) ^ 0. Further the dependence of the minimum value (denoted by |Im(a;o)|) of |Im(wq)| on 
the parameters can be summarized as follows. (1) The value of |Im(u;o)| decreases as Xc increases, and it becomes 
equal to ki for a certain value of Xc in the range 1 < Xc < 10 for the parameter values used in Fig. El (2) Even for 
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Xc = 1 it can become smaller than ki as the ratio Uoo/ni increases. For example, we obtain |Im(wo)| — 0.183 for 
Uoo = 40, Ki = 0.5 and a;c = ^ = 1. We can claim that the conditions Re(wq) ^ 0 and |Im(a;q)| < ki are satisfied in a 
wide range of the parameter values. 
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